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Abstract 

We study the asymptotic behavior, as p — > oo, of discrete measures 
on S*™ -1 that are induced by radially projecting point masses concen- 
trated on the integral lattice-points within dilates pD of a compact 
body D, for various classes of D. The results depend sensitively on 
the differential geometric properties of 3D. 

1 Introduction 

Recently, Douglas, Shiftman, and Zelditch, motivated by considerations aris- 
ing from the vacuum selection problem in string theory, have investigated 
questions in physics, aspects of which involve the equidistribution of radial 
projections of integral lattice-points on certain hypersurfaces in it" 1-1 , e.g., 
onto level surfaces of quadratic forms 2 . One question of this kind which 
they take up concerns the study of the asymptotics, as p — > oo, of the fam- 
ily of discrete measures on the surface d£ of an ellipsoid £ resulting from 
summing unit-weighted radial projections onto d£ of the lattice-points in 
p£ — {0}. If, for example, £ = S"™" 1 , this becomes: Does the action of 
the measure on S 11 " 1 which results from summing unit-weighted projections 
onto S n ~ l of lattice-points N € I/ 1 — {0}, with \N\ < p, and then rescaling 
the result by n/p n , tend, on smooth functions, to Lebesgue measure on 5 n_1 
as p — ► oo, and if so, how rapidly? This last question can be regarded as 
an instance of the general question of approximating smooth measures on 
5 n_1 by discrete measures, and in this paper, we will examine this question, 
which was suggested by questions arising in [2], which can be recast in this 
form. Depending on the measure being approximated, issues associated with 
the curvature of a specific surface associated with the measure can play a 
role, in particular the presence of zones of zero curvature on this associated 
surface, which we will illustrate by several examples. 



1 



It is a pleasure to mention that our interest in these questions arose from 
a conversation with Zelditch about his above-cited recent work with Douglas 
and Shiffman. 

We will begin the paper with a description of a general method, which 
can then be applied to various instances of the problem under consideration, 
in particular, to the zero curvature case, which is not discussed in [2J. An 
interested reader could probably glean the required background from a close 
reading of one or more of the papers [2], [7], ^U], as well as others, but in our 
view, expository clarity makes it very desirable to begin with such an exposi- 
tion, which we have written to lead as smoothly as possible into the notation 
and approaches of the various papers to which we refer. In particular, since 
detailed treatments of instances of the zero curvature cases are in general 
quite prolonged, and very closely mimic existing discussions in the litera- 
ture for the corresponding classical constant-density lattice-point problem, 
we have confined our discussion of these examples to general descriptions 
of how the overall techniques described in this paper can be adapted to 
closely follow treatments of the corresponding constant-density cases in the 
literature. 

2 

Suppose dp is a Borel measure on S n ~ x having a continuous, positive, piece- 
wise smooth density function m with respect to Lebesgue measure. If dp, as 
in the example above, happens to be Lebesgue measure, i.e., m = m(8) = 1, 
the corresponding family of discrete measures is parametrized by p > 0, and 
supported, as indicated above, on the unit-weighted radial projections of 
the non-zero integral lattice-points N, for which \N\ < p. Neglecting rescal- 
ing, the natural counterpart for a general measure on S"" -1 having positive 
density m with respect to Lebesgue measure is the family dT p of discrete 
measures supported on the unit-weighted radial projections of the non-zero 
integral lattice-points in pD, where D is the compact set whose boundary is 
given in polar coordinates by r = (m(#)) 1//n . We assume that m(6), or equiv- 
alently dD, is sufficiently regular so that the divergence theorem is valid for 
D - for example, dD could be smooth, or the boundary of a polyhedron. 

3 

We begin by noting that the effect of the above projection measure on a 
smooth function f{9) on S n ~ l is identical to that of the lattice-point count 
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over pD — {0}, weighted by the weight-zero homogeneous extension F of 
/ to R n — {0}. There are various analytical approaches to the estimation 
of such sums, e.g., Riesz means, as in [3], [3], etc., or convolution smooth- 
ing techniques, as in, for example, [2j and {?]• We will employ a convolution 
smoothing technique, since it is relatively simple to implement and describe. 
In our outline of this approach, we will temporarily assume, for minor tech- 
nical reasons, that D has been modified by the removal of a neighborhood of 
the origin. Then, assuming for the moment that / is positive, and denoting 
its homogeneous extension to R n by F, the above lattice-point count Nd(p) 
in pD can be overestimated by summing over the weighted lattice-points 
in a slight expansion E(p) of pD and underestimated by summing over the 
weighted lattice-points in a slight contraction C(p) of pD. The Poisson sum- 
mation formula provides a natural analytic approach for handling such sums, 
but since the function F is not smooth when multiplied by the indicator 
functions of E(p) or C(p), i.e., at the boundaries of E(p) and C(p), we will 
employ a slight refinement of the above idea, which can be approximately 
described by saying that we convolve F, restricted to E(p) and C(p), respec- 
tively, with a smooth compactly supported approximate delta function 8 e , 
whose support is chosen to lie within a ball of radius e(p) which is sufficiently 
small so that the convolutions are very close to the function F restricted to 
pD, except in small neighborhoods of the boundary. If the support of 5 t is 
of sufficiently small diameter, it will then be exactly true, if F is constant, 
and approximately true otherwise, that 5 e * i*b(p) — E p o < 5 e * F E t p \, where 
Fq denotes the product of F with the indicator function \& of S. If F is 
not constant, the deviation from correctness of the above inequality can be 
quantified in terms of a bound for the directional derivatives of F and the 
diameter of the support of 5. In view of this, if, on E(p) and C(p) respec- 
tively, we modify F E r p \ and ^b(p) by the addition and subtraction of a suit- 
ably small quantity whose size depends on p, the above inequality becomes 
correct, and can be exploited, following the method of [7j, which treated 
the case of constant density, to derive estimates for Nd (p) . If / is of mixed 
sign, we can express it as the difference of two positive smooth functions, 
proceed as above, and later combine the estimates for both. It is almost 
immediate that any smooth function is expressible as the difference of two 
positive smooth functions, whose size and the size of whose derivatives are 
comparable to those of the original function. For instance, we could define 
/ = /+-/-, where /+ = / + c, and /~ = c, where c > \mm 6eS n-if(6)\. 

We now show in greater detail how to carry out this program. The 
outcome in a general sense will be that after rescaling by n/p n , discrete 
sphere measures which arise in the above way will converge to dp on S 1 ™ -1 
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at the same rate as that at which l/p n times the standard unit-density 
lattice-point sum for pD converges to the measure of the set D. The specific 
sense in which this is the case will be developed more precisely below. 

4 

We will retain the notation of Section (J2J). As remarked there, the task is 
equivalent to the estimation, as p — > oo, of the lattice-point count over pD, 
with lattice-points counted with density F, where F is a smooth homoge- 
neous function of weight on R n — {0} (if F is constant, it is smooth at 
the origin as well). Unless F is constant, the origin causes a small, easily 
surmounted technical difficulty, which can be handled in various ways. We 
will, following Douglas, Shiffman, and Zelditch, deal with this by estimating 
the weighted lattice-point count in the shell pS = pD — ^D, which is a dilate 
of the basic shell S = D — \T). Since this estimate will be uniform in the 
dilation parameter, this automatically leads to estimates for the count for 
\pS, \pS, ^pS,. . . . We then obtain the desired result for pD by adding up 
the results for pS, \pS, jpS, ^pS, .... The weighted lattice-point count for 
pS will be estimated along the lines outlined above, i.e., by bracketing the 
true count between appropriate convolutions. 

5 

We now pass to the details. 

We will use the following lemma (cf. pp. 262-263), whose role is to 
show that the Fourier transform g(y) of a smooth function g{x) on a domain 
6 for which the divergence theorem is valid can be expressed as the Fourier 
transform of a smooth function on d&, and that the process of transfer to 
the boundary picks up a factor of l/\y\- 

Lemma 1. Suppose g(x) is a smooth function on R n , and y a fixed 
vector in R n . Then there exists a smooth vector field F(x) on R n , such that 

1. div [(2vri|y|)- 1 e 27ri ( :r '^F(x)] = g(x) 

2. The derivatives up to order k of the components ofF(x) can be bounded, 
independently of y, in terms of bounds for the derivatives up to order 
k + l of g{x). 

Proof. 
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This is proved in Lemma 3 of jSJ. In particular, since, setting x = 
(xi,...,x n ), y = (yi,...,y n ), and F = (F 1 ,...,F n ), the divergence of 
(2m\ y \)~ 1 e 27ri ^F(x) is 

e ^(».»)[(( !/1 /[ y |)F 1 + (2m\y\y 1 dF 1 /dx 1 ) + ... 
+ ((y n /\y\)F n + {2iri\y\)- l dF n /dx n )} , 

= e 2m ^[(y/\y\,F) + (27Ti\ y \)- 1 dWF], 
the requirements of the present lemma will be satisfied if 

1. divF = 

2. (y/\y\,F)=g 

with the derivatives of the Fj's bounded as indicated above. This is, however, 
precisely the assertion of Lemma 3 of jS], taking = y/\y\ in the notation 
of that lemma. □ 

Corollary. The first assertion of the lemma immediately implies, by 
the divergence theorem, that 

f e 2wiix ' y) g{x)dx = (27ri\y\)- 1 f e 2wi{x ' y) (F(x), n{x)) ds x , 
Je Jd& 

where n(x) is the exterior normal to d<5 at x. This gives the desired expres- 
sion of the Fourier transform over 6 in terms of a Fourier transform over 
d&. 

As previously indicated, we will assume that F > 0, since the general 
case can be deduced from this one by taking differences of the resulting 
estimates. For a small positive parameter e(p), which will depend on p, and 
which we will sometimes simply write as e for short, we define, respectively, 
a slight expansion E{p) of pS, and a slight contraction C{p) of pS, by setting 
E(p) = (p + e)D - (f - e)D, and C(p) = (p - e)D - (f + e)D. Since we 
are interested in the asymptotics of the weighted lattice-point count over 
pS as p — » oo, we may assume, by replacing S by a sufficiently large dilate 
of its original self if necessary, that for small e(p), E(p) contains an e(p)- 
neighborhood of pS, and that pS contains an e(/9)-neighborhood of C(p). 

Throughout the following, our basic approach is via the Poisson sum- 
mation formula, and is essentially the same technique that was employed 
in [7], the only difference being that now the weight function is no longer a 
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constant. Since N$(p) = ^2ngpS a purely formal application of the 

Poisson summation formula leads to 



(1) 



as an analytic expression for the weighted lattice-point count, where i*e 
denotes the Fourier transform of Fq. 

As previously mentioned, there are convergence difficulties, since the 
function which is equal to F(x) on S and equal to zero on the complement 
of S is discontinuous at the boundary of S unless F = 0. This problem is of 
course present in the constant density case as well, and it can be addressed, 
as was done in [7] for that case, and as we have previously indicated, by 
slightly expanding and contracting pS, into E{p) and C(p), respectively, 
and then smoothing the indicator functions of E(p) and C{p) by convolving 
them with a smooth approximate delta function 5 e (x), whose support lies 
in a ball of diameter e(p). Then if F is constant, 



and since the functions being summed are smooth with compact support, 
the Poisson summation formula can be applied to both sides of the above 
inequality to obtain asymptotics for N$(p) (cf. [Jj). 

In the case of non-constant F(x), a minor technical issue connected with 
this approach arises from the fact that F{x) is no longer necessarily equal 
to its convolution by 5 e (x). However, for large p, (5 e * F){x) is very close to 
F on a dilated shell p<5, since any first derivative of F is homogeneous of 
weight —1, and therefore the oscillation of F over the intersection of a ball 
of diameter e with p6 is at most of the order of (l/p)e(p). I.e., F is close 
to being constant on p6 as p becomes large, so, as previously indicated, the 
addition and subtraction of a suitable small quantity r\ p in inequality (J2J) for 
the constant density case results in a correct inequality for the general case: 



^ [6 e *(F c -TjpXc)KN)<Ns(p)< Yl [Se*(FE + V P XE)](N), 




(2) 



Nec(p) 



NeE(p) 



N<=C(p) 



N£E(p) 



or 
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(5 e *F c )(N)- (Se*V P Xc)(N)<N s (p) 
NeC(p) Nec(p) 

< J] (6 e *F E )(N) + (Se*V P XE)(N), (3) 

NeE(p) NeE(p) 

where x& denotes the indicator function of 6, and for notational simplicity 
we have written C and E for C(p) and E(p), respectively. 

Since r] p will be chosen to be of the order of (l/p)e(p), the analysis of 
the contribution to the above inequality that is attributable to the presence 
of the rj p terms on either side will be easily seen to lead to the presence, on 
either side, of a term of the order of p n ~ 1 e(p). The choice of e as a function 
of p will depend on the shape of D, and for the time being, we will simply 
carry it along as an unspecified function of p. 

When we come to the analysis of the sums in © , the heart of the matter 
lies in the asymptotics of the Fourier transform of F over C(p) and over E(p), 
which reduces, by an obvious rescaling, to the analysis of the asymptotics 
of the Fourier transform Fq of F over shells S which are very close to S for 
large p. 

Accordingly, we next take up the asymptotic analysis of 

F 6 (y)= [ F(x)e 2 ™^dx, 
J& 

where S = &(p) can stand for a rescaling by 1/p of either C{p) or E(p). 
By Lemma 1, 

[ e 2 *< x ^F{x)dx=(2m\y\)- x ! e 2wi(x ' y) {F(x),n{x)) ds x , (4) 
J& Jae 

where the size of the derivatives of the vector field F are controlled by the 
size of the derivatives of the homogeneous function F, which, since S is a 
shell, we can smooth at the origin without affecting its values in S. 

The boundary of & splits into two components, which are analyzed sim- 
ilarly. Thus, the asymptotic analysis of 

/ F{x)e 2 ^ y) dx 
J& 

amounts to the study of the asymptotic behavior of integrals 

J g{x)e 2 ^ y) dx, 
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where <S is a hypersurface in R n of the considered type, and g is a smooth 
function on S. This is a highly studied subject (cf. [3], [7j, [S], et al.). The 
results are very dependent on the geometry of S, although the case in which 
S has everywhere positive Gaussian curvature is easy to analyze, using, for 
example, a straightforward application of stationary phase, or even a simple 
calculus argument. It is well-known that in this case (cf. [3]), 

j g{x)e 2 ^ y) dx < \y\~ {n - 1)/2 , (5) 

where the implied constant depends on bounds for the derivatives of g, 
which, for the purposes of this discussion, can be thought of as a function 
defined in a neighborhood of S. By the corollary to Lemma 1, it follows 
that in the positive curvature case, for an n-dimensional shell ©, 

[ e w ^F(i) dx < [yr (n+1)/2 . (6) 
Je 

Although the positive curvature case has already been discussed in [2J, we 
will use this case to illustrate our approach, since it provides an exceptionally 
simple paradigmatic description of techniques of considerably more general 
applicability. We will then describe modifications that are required in illus- 
trative cases for which dD does not have everywhere positive curvature. In 
general, specific applications require individual adaptations of the method, 
but these adaptations can often be easily inferred from existing treatments 
of the constant-density case. 

The technique closely follows that of jjj, and so will be simply outlined 
here. Denoting as above by Ns{p) the weighted lattice-point count in pS, 
we will indicate how to use the Fourier transform bound to estimate the 
right side of the inequality (jSJ). The left side is handled in a similar way. 
The Poisson summation formula can be applied to the right side of ©• We 
first take up the term 

Y (5 e *F E )(N), 

N&E{p) 

which, by the Poisson summation formula equals 

Y Un)f e{p) (n), 

N€E(p) 

= <5(0)F(0) + Y! UN)F E{p) (N), 
NeE(p) 
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where the prime on the summation sign means that the origin is omitted 
from the sum. 

Since 8(0) = 1, the last quantity can be rewritten as 

/ F(x)dx+ Y* 8 € (N) F E(p) (N) . 

JE(p) 



NeE(p) 



If in the above we replace the term f E r p \ F{x) dx by J pS F(x) dx, the 
error corresponding to the replacement can be estimated by noting that 



/ F(x) dx - / F(x)dx = I 
JE(p) JpS Je 



F(x) dx 



E(p)-pS 

« vol(E(p)-pS) 



n-l 



« ep 

The remaining term on the right side of ©, namely 

£ (5 e * VpX E)(N), 

NtE(p) 

is clearly <C (rj p times the volume of E p ), i.e., of order rj p p r ' 
We can thus rewrite the right side of (J21) as 



/ 

JpS 



F(x) dx + 0(ep n 



~ 1 +rip n \ 



+ Un)f e{p) (n) 

NeE(p) 



(7) 



The last sum is handled in exactly the same way as in the well-known 
case in which F is constant. Namely, the necessary estimate on F E ( P ){N) 
is provided by ©, and is identical with the estimate for the constant case. 
As in that case, integration by parts shows that the term 5 e (N) is 1/(1 + 
e|iV|) fc , for any fixed integer k. The sum is estimated by splitting it into 
two parts, over N for which |7V| < 1/e, and over N for which \N\ > 1/e, 
respectively. Each of these parts is estimated by comparison with an integral, 
using a sufficiently high value of k to produce convergence, and the estimate 
(jHJ) for the Fourier transform. The result is then minimized by choosing e 



to balance the estimates, which leads to the choice of e 



P 



>-l)/(n+l) 



which in turn results in an estimate of p( n 1 )( n /( n + 1 )) for the infinite sum 
in Q. The above choice of e results in the same estimate for the term of 
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order ep n ~ l in (|7|). and since rj is of the order of p _1 e, we obtain the same 
estimate for the term np n in (|7|). 

Since a similar analysis can be applied to the left side of the inequality 
(j3J), this shows that in the positive curvature case, the weighted lattice-point 
count over pS equals 

/ F(x) dx + O^- 1 ) ) . 

J P s 

In order to analyze the weighted lattice-point count over pD, we add up 
repetitions of the above quantity, corresponding to p, p/2, p/2 2 , . . . , p/2 k , 
where k is taken to be of the order of log 2 p, so that p/2 k < 1. This takes 
account of all lattice-points in pD except perhaps for a fixed finite number 
near the origin, and since these latter have no effect on the asymptotics, we 
find that 



N D (p) = I F{x)dx + Y J 0{{p/2 j )^ l ^ n ^ n+1 ^) 

J P D j=Q 

= p n I F{x)dx + yO((p/2 j )^ n - 1 ^ n ^ n+1 ^). 



j=0 

Since the implied constants in the O terms are uniformly bounded, this 
implies that 

N D {p) = p n [ F(x) dx + 0( /9 («- 1 )(™/(«+ 1 ))) . 

JD 

Now it is a consequence of equation (4) of [3] that 



F(x)dx = (l/n) f(0)m(9)d0, (8) 

D Js™- 1 

so we obtain the following equivalent form of the result of Douglas, Shiffman, 
and Zelditch in the positive curvature case: 

Theorem, (cf. [2]) With notation as above, 

[ f(6)m(6) d9 = (n/p n )N D (p) + 0{p- 2n /^) . 
JS"- 1 
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This is an estimate for the rapidity of the convergence of the discrete mea- 
sures (n/p n )dT p to dfi = m{6)d6. The estimate, which was given in a 
somewhat different form in coincides with the usual scaled error term 
for the standard lattice-point problem for bodies with boundaries having 
strictly positive Gaussian curvature (cf. [3]). 

To justify we note that from (4) of [3] it follows that 

F{x)dx= f t n ~ l dt I F(tx)(x,n(x))ds x , 

D Jo JdD 

where n(x) is the outward normal to dD. Since F is homogeneous of weight 
0, the double integral equals 



(1/n) / F(x) (x,n(x))ds x , 

JdD 

and coordinatizing dD by 5 n_1 via the radial map, this becomes 

/(0)*(0) (x,n(x))d9, 

S n-1 

where is the Radon-Nikodym derivative ds x /d6, and x is regarded as a 
function of via the radial map. It follows easily from elementary geometric 
considerations that 

$(0) = \ x \ n /(x,n{x)), 



so 



F(x)(x,n(x))ds x = / f(9)$(9)(x,n(x))d6 

dD Js 71 - 1 



ii.-i 



\x\ n f{9)d0. 

But on S n ~ l , \x\ = (m(#)) 1//n , so the last integral equals 

f(9)m(0)d9, 



n-l 



from which (03) follows. 
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6 The Case in which Curvature can Vanish 



Classical lattice-point asymptotics for dilates pD of a body D for which the 
curvature of dD is not always positive can be quite intricate, and are highly 
dependent on the manner in which the curvature of dD vanishes, as well 
as on the placement of D in relation to the integer lattice. For example, 
in R 2 , if the curvature of dD vanishes to finite order at a finite number of 
points, the classical lattice-point asymptotics depend on the order to which 
the curvature vanishes at the points in question, as well as on Diophantine 
properties of the normal vectors to dD at those points (cf. [7j). The situation 
can become much more complicated in higher dimensions, and nuances of 
this type similarly affect convergence rates of the discrete measures with 
which we are dealing in this paper. Since, once one has the necessary Fourier 
transform asymptotics, methods for dealing with these issues closely mimic 
those for the classical case, we will content ourselves with briefly indicating 
what happens in a few interesting representative cases. 

As indicated above, the central analytic issue is the detailed asymptotics 
of the Fourier transform of smooth functions on dD. In the case of every- 
where positive curvature, one has the previously mentioned result that 

/ g(x)e 2m{x ^ dx « \y\- {n ~ l)/2 , 
Jav 

where the estimate does not depend on the directional component of y. 

This estimate is generally false if the curvature of dD vanishes on some 
non-void subset of dD, and a useful description of what happens can be 
quite complicated, for example, if 

dD = {(x u . . . , x n ) | x x 2k + ■■■ + x n 2k = 1} , 

or if D is a polyhedron, (cf. [Bj and 
If we write 

f e ^(x,y) dx 
Jd 

in polar coordinates as \I/(r, </>), one quite general fact along these lines is 
given by Theorem 1 of [Bi. Namely, if dD is real-analytic and D is convex, 
then the function 

A(4>) = supr( n+1 )/ 2 ^(r,</>) 

r 

is in I/^S" -1 ), for some p > 2. Thus, under conditions of considerable gen- 
erality, the Fourier transform asymptotics coincide with those of the positive 
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curvature case, up to multiplication by a function in L^S 1 ™ -1 ). The con- 
vexity hypothesis is unnecessary in R 2 , and possibly in higher dimensions, 
although this is not generally known. The requirement of real-analyticity 
can be replaced by somewhat weaker hypotheses (cf. m]). For later related 
papers, cf. PQ, [To] , 

Example 1 

The arguments by which the above-mentioned theorem is proved can be 
uneventfully applied to obtain the following counterpart of (JSJ), in the case, 
for example, in which D is convex and dD is real-analytic (and under weaker 
hypotheses in R 2 ): 

[ g{x)e 2wi( - x ^ dx <^ A(</>)br (n ~ 1)/2 , (9) 
J 3D 

where is in L^S"™ -1 ), for some p > 2. This implies, by a straight- 

forward application of the techniques of this paper, that if we modify the 
definition of the discrete measures T p by replacing Z n by its image under 
the action of an element 7 of SO(n), call the resulting measures ^(7), and 
denote the sum corresponding to Nr,{p) by N£>{p, 7), that we obtain the 
following result: 

Theorem A (vanishing CURVATURE case). With notation as above, 
[ |i? D (p, 7 )|rf7«0(p- 2n/(n+1) ), 

iSO(n) 

where 

R D (p,l)= I f(e)m(e)d6-(n/p n )N D (p n ). 
JS"- 1 

I.e., even for a quite general version of the zero curvature case, the error 
estimate for the positive curvature case holds for the L 1 norm over the 
rotation group of the errors for the "rotated" measures. There are, of course, 
various consequences corresponding to other LP norms as well. 

Example 2 

Our next example is the special case in which 

dD = {(xi, . . . , x n ) I Xl 2k + ■■■ + x 2k = 1} . 
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The relevant estimate on the Fourier transform is given by Theorem 2 of [S] , 
which states that for sufficiently smooth g, if we express 



Jdv 

in polar coordinates as 9(r, (f>) (cj) € S n , 4> = (<f>\, . . . , then on the 
set of points (r, c/>) for which exactly j of the (f)*'s vanish, 



l)/(2k — 1), and ctj = (j/2k) + (n—j — l)/2. This is, of course, in some sense a 
special case of ©, but it gives considerably more detailed information about 
the asymptotics of the Fourier transform for this case. 

The classical lattice-point problem for this case is discussed in 6 . The 
principal result is that the error term is of order p R , where R = max(>l, B) 
with A = (2k — l)(n — l)/2k and B = n(n — l)/(n + l). The estimate is best 
possible if A > B. This result is obtained by separately analyzing groups 
of lattice-points on the Fourier transform side of the Poisson summation 
formula, where the grouping is arranged into sets defined by the vanishing 
of a specific number of the </>*'s. That the result is sometimes optimal is a 
consequence of the observation that if A > B, the contribution to the Fourier 
transform side of the Poisson summation formula coming from lattice-points 
on axes defined by the vanishing of all but one of the (/>*'s constitutes the 
major contribution to the error, and that the behavior of this contribution 
can be analyzed by a stationary phase argument. 

The adaptation of this result to the present context is straightforward, 
and the arguments are nearly identical to those for the original result, so 
as before, we will content ourselves with a statement of the theorem in the 
present context, using notation as above. The result in the present context 
is: 

Theorem B (zero curvature case). With notation as above, if 





(10) 



dD = {( Xl ,...,x n ) \x 1 2k + --- + x. 



2k 



1} 



n 



then 




and if A > B, this is best possible. 
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Example 3 



As a final example, we will consider the polyhedral case, and in the interests 
of expository and combinatorial simplicity, describe a typical result for the 
2-dimensional case. 

In the classical constant-density lattice-point problem, if at least one 
of the perpendicular vectors to a face of a compact polyhedron has rational 
coordinates, there are an infinite number of pi — > oo from which an infinitesi- 
mal displacement results in a modification of the lattice-point count of order 
p n ~ l , so in this circumstance the error estimate is of true order p n ~ x . Since 
a simple estimate of Gauss shows that the error term is always <C p n_1 , 
polyhedra can be worst possible cases for lattice-point error asymptotics. 

Paradoxically, this situation is not generic for polyhedra, as was noticed 
by Khintchine 4 in the 2-dimensional case. Khintchine's result is that the 
error estimate corresponding to almost any rotation of the integer lattice Z 2 
is in fact extremely small. Specifically, for any e > 0, it is almost always 
< log 1+e p- Later work [7j, 0, QZ3, EB, EH, PH has considered 

various aspects of the n-dimensional well as additional features and 

refinements of the 2-dimensional case. 

As in the previous examples, methods for the classical constant-density 
lattice-point problem carry over with very little change to our context. As 
a typical illustrative representative of what can be expected, we begin by 
recalling the result, due to Skigranov [TT] . that if D is an algebraic polygon, 
(one for which the ratios of the direction numbers of the normals to its 
faces are all algebraic of degree > 2), then the classical lattice-point error 
term is <C p e , for any e > 0. This result is also a special case of theorems 
in later papers E3 an d |12j . It can be obtained by using a more detailed 
form of the estimate © , as was done in the previous example (|l(Jj) , but one 
which is adapted to the particular case in which D is a polygon (cf. [5]). 
With such an estimate in hand, the lattice-points on the Fourier transform 
side of the Poisson summation formula are split into two groups: those in 
finite- width bands surrounding the normal vectors to the sides of D, and all 
the rest. The contribution from the lattice-points exterior to the bands can 
be estimated by comparison with an integral, while the series arising from 
the contributions from lattice-points within bands is estimated by using 
Diophantine properties of the ratios of direction numbers associated to the 
corresponding normals. Since the relevant estimate for the Fourier transform 
of D is singular at these directions, the poor approximability of these ratios, 
which is a consequence of Roth's Theorem, is crucial (cf. e.g., p. 858 for 
a similar argument). The corresponding result, in the context of the present 
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paper and expressed in the notation of this paper, is that for a polygon of 
the above type, 

Theorem C (polygonal case). With notation as above, 



7 Conclusion 

We have described a general method for describing the accuracy with which a 
large class of measures on S n can be approximated by a naturally associated 
family of discrete measures. The case in which dD, in the notation of this 
paper, has everywhere positive curvature has been previously studied in [2], 
and is taken in this paper as a basic template for the description of a general 
approach to such problems, in particular, cases involving zero curvature. 
As in the classical constant-density lattice-point problem, there are special 
instances of the positive curvature case, e.g., arithmetically defined positive 
definite quadratic forms, in which the general estimate can be improved by 
exploiting the underlying arithmetic character of the associated surface. In 
the case in which dD contains subsets on which the curvature vanishes, the 
situation becomes vastly more intricate, although there are general results, 
e.g., along the lines of the first of our three examples in the zero curvature 
case. There are also arithmetic instances of the zero curvature case having 
special features, as in the second of our three examples. One can give general 
results for polyhedra as well, e.g., along the lines of 02] ^E]- The 
methods will in general mimic those for the constant-density case, once 
one is in possession of the appropriate Fourier transform asymptotics (for 
the asymptotics in the polyhedral case, cf. [!) ). A kind of meta-conclusion 
is that in general, the derivable asymptotics associated with the presently 
considered class of problems coincide with the corresponding results for the 
classical constant-density case, by virtue of the fact that the relevant Fourier 
transform asymptotics are effectively identical. 
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